High-Order Variational Calculation for the Frequency of Time-Periodic Solutions 
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We develop a convergent variational perturbation theory for the frequency of time-periodic solu- 
tions of nonlinear dynamical systems. The power of the theory is illustrated by applying it to the 
■ Duffing oscillator. 
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which is also known as the Duffing equation ||. Here the dot abbreviates the derivative with respect to the time t, 
lu denotes the harmonic frequency, and g > stands for the coupling constant. In the following we solve (|l]) for the 
5^ initial values 

x(0) = l, 4(0) =0, (2) 
and determine the frequency u of the resulting periodic motion by using variational perturbation theory. In Section || 
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O . I. INTRODUCTION 

o 

Perturbative treatments of physical problems provide us with divergent power series in some coupling constant g. 
Typically, the perturbation coefficients grow factorially, so that they have a zero radius of convergence. If the signs of 
the perturbation coefficients alternate, there exist various resummation schemes which help us to obtain finite results 
for all values of the coupling constant g, even in the strong-coupling limit g — > 00 (for an overview see Chap. 16 
of Ref. Q). Most successful is variational perturbation theory which was recently developed (2[ ||| as a systematic 
extension of the variational approach of Feynman and Kleinert |ij . Initially, this theory was set up for calculating 
the effective classical potential in quantum statistics. It has been thoroughly tested for the ground-state energy of 
the anharmonic oscillator and shown to converge exponentially fast and uniformly to the correct result |^|, [|. This 
£S) I success has led to applications to divergent series in other branches of theoretical physics g]. Most spectacular was 
the success in calculating the most accurate critical exponents of the 4 -theory without using re-normalization group 
1 methods fj], fj}. 

Oh' 

In this paper we extend variational perturbation theory by developing an exponentially fast converging variational 
perturbation theory for the frequency of time-periodic solutions of nonlinear dynamical systems. As a simple but 
■ nontrivial model we consider the one-dimensional anharmonic oscillator with the equation of motion 



x(t)+Lulx{t)+gx 3 (t) = 0, (1) 



we calculate the frequency w as a power series of the coupling constant g. Section III then elaborates the variational 
resummation of this weak-coupling series so that the frequency lu can be determined for all values of the coupling 
constant g including the strong-coupling limit g — > 00. 

II. PERTURBATION THEORY 

We start by solving the initial value problem ([!]) and (Q) perturbatively to high orders. 

A. Poincare-Lindstedt Method 



We assume for a sufficiently small coupling constant g that the solution x(t) has the asymptotic representation 

x(t) =x (t) + Xl (t)g + ... . (3) 
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A systematic standard procedure to obtain such an asymptotic series for a periodic solution 

x(t) = x (t + —j 



(4) 



is provided by the Poincare-Lindstedt method ||, [l0|]. There one explicitly takes into account that the unperturbed 
frequency loq is shifted to the frequency wbya nonzero coupling constant g. One performs a rescaling of time according 
to 



i(0 =*[-. 



and introduces the new variable 



This converts the periodicity condition (|J) to 

<?(£) = 9 (£ + 2tt) 
and transforms the original initial value problem (^) and (||) to 

u> 2 q"(0+"U(0+9Q 3 (0 = 0, 9(0) = 1, </(0)=Q, 



(5) 



(6) 



(7) 
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where the prime indicates the derivative with respect to the dimensionless new time variable £ . Since the coupling 
constant g is supposed to be small, we can expand the frequency lj and the period solution in powers of g 
according to 



= E w « w °(5) ' 



71=0 



(9) 
(10) 



Due to this ansatz the expansion coefficients w n and g„(£) are dimensionless. Inserting (|^) and (|T^) in the initial value 
problem (||) and comparing equal powers in the coupling constant g leads for n = 1,2,... to the following recursive 
set of ordinary differential equations: 



q'n(0 + g»(0 = /»(£) , g«(o) = ^(o) = o , 

where the inhomogeneity / n (£) is given by 



(11) 



n— 1 n—m 



n—1 n—m — 1 



/ n (0 = -2«, n 5o '(0 - 2 £ - X) E ^ C TO -«(0 -EE «(0 ?»-™-i-i(0 • (12) 



7 = 1 



m=l 1=1 



m=0 Z=0 



This is solved starting from 



wo = 1 , go(0 = cos £- 



(13) 



In the nth integration process we proceed in three steps. At first, we calculate the inhomogeneity / n (£) according to 
( |l2| ) and expand it in a Fourier series which turns out to be of the following form: 



/»(£) J^fn.k cos(2£: + l)£. 



(14) 



fe=0 



Second, we prevent a secular term in q n (£) from solving (gdj) by demanding the condition 

/n,l = 0, 



(15) 



from which the expansion coefficient w n is uniquely determined. Third, the initial value problem (11) is solved by 

J 7i , k 



E {2 k + 1)2-1 



fc=i 



fn,k 

(2k + iy 



cos(2fc+ 1)£. 



(16) 
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n 




1 


3 
8 


11 


35H276321347 
562949953421312 


2 


21 


12 


401225915283063 




~256 


72057594037927936 


3 


81 


13 


2892201453147555 


2048 


576460752303423488 


4 


6549 


14 


84053106665670789 




262144 


18446744073709551616 


5 


37737 


15 


614845335384090729 


2097152 


147573952589676412928 


6 


936183 


16 


1158192705499996341141 


67108864 


302231454903657293676544 


7 


6077907 


17 


8566538482894401288225 


536870912 


2417851639229258349412352 


8 


2604833685 


18 


254612814518190043882263 




274877906944 


77371252455336267181195264 


9 


17839453041 


19 


18996276910402y2362960331 


2199023255552 


618970019642690137449562112 


10 


497158650207 


20 


227596989316436230247319519 


70368744177664 


79228162514264337593543950336 



TABLE I: The first 20 dimensionless weak-coupling coefficients w n for the frequency to of the Duffing oscillator. 



Using a computer algebra program we obtain in this way the perturbation expansions for both the frequency 



and the periodic solution 



u = lu + —g - T^r^a +■■■ ( 17 



1 1 



x(t) = cos ut + ( — or| 3 cosa;t + 3 cos3a;£ 



23 3 

j cos wt t cos 3wi H t cos 5wi I q z + . . . . (18) 

4 i no, .4 i no/, ,4 ( » 1 v / 



1024^ 128^ 1024w 



Tab. | shows the first 20 weak-coupling coefficients w n of the frequency w. 



B. Analytical Expression for the Frequency 

Remarkably, the frequency u> of the Duffing oscillator (^) can be determined exactly as a function of the coupling 
constant g. To this end we multiply ([j]) by x(t) and integrate once, taking into account the initial values (^): 

h 2 {t) + X -loI x 2 (t) + \gx\t) = 1 u; 2 + i «? , (19) 

Separating the variables in the energy conservation (|l^) and integrating over a quarter of the oscillator period, we get 

1 ' (20) 
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The integral can be explicitly performed by using Eq. (3.152.4) in Ref. [Tl| 



2F 



(21) 



,2'V 2 (C0+S) 

where F denotes the elliptic integral of the first kind which is defined in Eq. 



.111.2) of Ref. p] 



da 



V 7 ! - k 2 



(22) 



sin a 



In the weak-coupling limit g — > 0, we recover from ( |2l| ) the perturbation series (^) of the previous section by taking 
into account (8.113.1) of Ref. 11 . However, the exact result (21 ) also allows to find the strong-coupling limit g — > oo: 



fe + b 



5 



(.9 oo) • 



The leading strong-coupling coefficient 6o has the numerical value 

7T 



6o = 



2.F 1 



7T 1 

2'7! 



= 0.8472130847939790866 .... 



(23) 



(24) 



In Figure [l] we compare the full function (pTh with the successive divergent weak-coupling expansions (|^) and with 
the convergent strong-coupling expansions (p3) . We observe that the full function represents the envelope to all weak- 
and strong-coupling expansions. 



III. VARIATIONAL PERTURBATION THEORY 



Let us now see how well we can reproduce the strong-coupling result by resumming the weak-coupling series @ of 
the frequency to with the help of variational perturbation theory ^| . In this way we are able obtain approximations 
of the frequency lu in good agreement with the exact result for all values of the coupling constant g. In particular, 
we carry the strong-coupling limit g — > oo of the frequency lu to high orders in order to investigate in detail the 
convergence of the variational results. 



A. General Procedure 



We start with the weak-coupling expansion (ft) truncated at order N: 



N 



,(N) 



2n g-n 



The we introduce the variational parameter Q by Kleinert's square root trick 

lu = J Q 2 +iol~ Q 2 = Qy/l+gr, 

where the abbreviation r is defined by 



gW 



(25) 



(26) 



(27) 



and we ignore for the moment that r depends on g and regard it as a constant. Substituting (|26|) into the truncated 
weak-coupling series (|25|), we obtain 



N 



(g,n) = Y,^ 1 ~ 2n (l + gr) 1/2 - n g n 



(28) 



5 




lg 9 



FIG. 1: Logarithmic plot of the full function (211) (solid curve) versus the coupling constant g compared with the truncated 
successive divergent perturbation expansions ( |) (dashed curves, the corresponding orders are labeled with 1 ... 9) and the 
partial sums of the convergent strong-coupling expansions ( p3| ) (dotted curves, the corresponding orders are labeled with 
1...9). 



The factor (1 + gr) a with a = 1/2 — n is then expanded up to the order N — n 

N-n / \ 

a 



k=0 



where the binomial coefficient is defined by the Gamma function T: 

T{a + 1) 



a 
k 



r(fc + i)r(a + fe + i) ' 

Thus the sum (|2§| ) is re-expanded including all powers of g up to the order g N 

N 

At the end we have inserted (p7j). 



N—n 

E 

k=0 



1/2- n 
k 



n 2 



(29) 



(30) 



(31) 



If we could consider the limit N — > oo in ( ^l| ) , the dependence on the artificially introduced variational parameter 
fi would drop out. Due to the truncation at the finite order TV, however, we obtain an explicit dependence on the 
variational parameter f2 in (|3l]) . This suggests to fix the yet undetermined variational parameter Q according to the 
principle of minimal sensitivity pM- Thus we try to find at first an extremum: 



an 



0. 



n=n<™)( 3 ) 

If this equation has no real solution (g) , then we look for a saddle point instead 

d 2 uj<- N \g,n) 



an 2 



o. 



(32) 



(33) 
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or more generally, for a real zero Q( N \g) of the lowest derivative with respect to the variational parameter. This 
optimal value ^ N \g) then leads via 

w W( ff)=w W( 5 ,nW( 5 )) (34) 

to an approximation of the frequency uj which turns out to lead to good results for all values of the coupling constant 
g. To lowest order, the optimal solution is unique. If there are several optimal solutions, we always choose the one 
which is closest to the optimal solution of the previous order. 



B. Lowest Orders 



Let us consider the lowest order approximations of variational perturbation theory explicitly. Truncating the weak 
coupling series (17) at the order N = 1 leads to 

3 



8uj 



Applying the general procedure as described in detail in the previous subsection, we obtain 



w «Q ? ,n) = ^ + 



n 



3.9 \ 1 

n 



From the condition (B3) we determine the variational parameter as 



(35) 



(36) 



(37) 



whose substitution into (pq) leads to the first-order approximation 



(38) 



In a similar way we proceed for the second order, where the truncated weak-coupling series ( |17j ) reads 

3 21 



CJ (2) = UJ + 



8U! 



9 



256wg 
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Applying the square root trick leads here to 



,(2), 



8' 







H 




16 J 





16 



1 



256 J n 3 ' 



(39) 



(40) 



It turns out that a/ 2 ) [g, f2) has no real extremum with respect to f2. Thus we have to look for a turning point instead. 
From condition (p3|) we obtain 



rt 2 \g)=u; 



\ 



35 , 21 5 2 



2w 2 32^4 



1 + 



(41) 



Inserting ( p~Ij ) in (|4^) leads to the second-order approximation 



uj {2) (g) = uj ■ 



1 , 3g | 153g 2 
•Jl 2560Jo 



Au 2 



35 , 21 3 2 



2uj\ 32u;4 



(42) 



In Fig. H we compare the first- and second-order variational approximation ( |38| ) and ( 42 ) for the frequency ui with 
the exact result (^l|). Notably, the first-order variational result is very good for all values of the coupling constant g, 
and the second-order leads to a substantial improvement of the accuracy. 
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C. Strong-Coupling Limit 

In order to quantify the accuracy of the variational approximations (pq ) and (0), we study now in particular 
the strong-coupling regime g — + oo. In first and second order we reproduce the expansion ( J23|) with the leading 
strong-coupling coefficient 

&« = ^ ra 0.866025 , 4 2) = ra 0.851895 . (43) 

o 2 224 

Comparing (|24|) with ([l3]) we see that the first and the second order of variational perturbation theory yields the 
strong-coupling coefficient ( p4[ ) within the accuracy of 2.2 % and 0.55 %, respectively. 

In order to obtain higher-order variational results for this strong-coupling coefficient (p4]) ; we proceed as follows. 
From the first- and second-order approximations (p7|) and ( ffl] ) we see that the variational parameter has a strong- 
coupling expansion of the form 

^ (iV) (fl)=^^ iV) + ^ + ^ + -J ■ (44) 
Inserting (Q) in (|3l]), we obtain the iVth order approximation for (H|) with 



N N-n / \ 

f(if)=EE 1/2 ; n (-in K ] ) 

n=0 fc=0 V fe / 



l-2n 



(45) 



The inner sum can be done by using Eq. (0.151.4) in Ref. |TT[: 



w 3 ) - ec-d*- f _ ^ /2 _; n ) «■* K )" 2 " ■ (46) 

n=0 \ / 

In order to optimize the variational parameter f2g^ we look again for an extremum 

O0 ° ^ I = (47) 

or for a saddle point 



d 2 b (N) (n {N} ) 

d b ° > =0. (48) 
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N 




N 


b (N) 


1 


0.86602540378443864676 


11 


0.84721311260106078088 


2 


0.85189520859585272618 


12 


0.84721309038427087031 


3 


0.84798320787226284162 


13 


. 8472 1 308733437656 102 


4 


0.84736735286736694523 


14 


0.84721308530703137833 


5 


0.84726277296604748829 


15 


0.84721308503241446175 


6 


0.84722291812428697005 


16 


0.84721308484231654612 


7 


0.84721687569394258505 


17 


0.84721308481682089873 


8 


0.84721383828896139276 


18 


0.84721308479862454760 


9 


0.84721340071349571092 


19 


0.84721308479620273029 


10 


0.84721314796371865932 


20 


. 8472 1 308479443254 1 39 



TABLE II: The first 20 variational results for the strong-coupling coefficient (|24| 



It turns out that an extremum exists for odd orders N, whereas even orders N lead to a saddle point. Tab. [n] shows 
the first 20 variational results for the strong-coupling coefficient (|4|). Note that the approximant fog 20 -* coincides 

already in 11 digits with The points of Figure 

linearly on N up to the order N = 100 according to 

\b N) 



show that the logarithm of the error \b N ^ — 6o|/^o depends 



bo\ 



-a-PN 



b 



(49) 



where the fit of the last 10 points leads to the quantities a — —6.7671 and (3 = —1.1113. Thus we have demonstrated 
that the variational approximations for the frequency of the Duffing oscillator converge exponentially fast. Note that 
the speed of convergence is considerably faster than the exponential convergence of the variational results for the 
ground-state energy of the anharmonic oscillator ||, [y] . 

IV. CONCLUSION 

We have demonstrated by the example of the Duffing oscillator how variational perturbation theory is successively 
applied to determine the frequency of time-periodic solutions of nonlinear dynamical systems. It remains to proceed 
along similar lines to treat also nonconservative systems with limit cycles like the van der Pol equation H. 
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FIG. 3: The points show the logarithmic plot of the error I&q^ — bo\/bo against the order JV, and the solid line represents a fit 
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